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Nonnegative Trigonometric Polynomials and 
Certain Rational Characteristic Functions 
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Let 0<6i<6 2 < • • • "On be n integers and 0<d 1 <Cd 2 <i . . . <^d m be m(m^ri) real 
numbers (not necessarily integers). Denote by g{6) the Vandermonde determinant formed 

m 
from the b\, b\, . . ., bl with the first row replaced by i — X, cos 6,0, where X,= II (1 — 6?/c© 

Jk=l 
for j=l, 2, . . ., n. The question whether g(0) is nonnegative for all values of 6 is closely 
connected with the problem whether certain rational functions are characteristic functions. 
Four configurations of the b\, . . ., b n , and d lf . . ., d m are studied, which lead to non- 
negative trigonometric polynomials. 

1. Introduction 

In this paper we study several nonnegative trigonometric polynomials that are connected 
with the question whether certain rational functions are Fourier transforms of distribution 
functions (characteristic functions) . 

We consider in the following a rational function <t>(t) w^hose poles and zeros have the same 
imaginary part a and assume further that all poles and all zeros have the order one. Let 
—ia, ivj, —ivj(j=l,2, . . . ,n) be the poles and —iw k , —iw k (k=l,2, . . . ,ra), m<n, be the 
zeros of 4>(t), where v j =a+ib j) w k =a+id k , and a>0, 0<6i<6 2 < • . • <ib n , 0<di<d 2 < . . . <id m . 
Such a function <j>(t) fulfills the necessary conditions that a rational function must satisfy in 
order to be the Fourier transform of a distribution function. These conditions were given 
in two earlier papers [1, 2]. 3 

In the following we denote by 



V k = 



b\\ 



h 2k 



bl 



&?'"- 



&!*" 



(1) 



the generalized Vandermonde determinant formed from the numbers b\, . . . , bl with the 

exponents k, 1, 2, . . . , (n— 1). We write further Aj for the minor of the element in the first 

n 
row and jth column of this determinant, so that Vj e =y y ,(—\V~ l b^ k A ,-. 

;'=i 

In order to decide whether the function <£(£) is a characteristic function, we decompose it 
into partial fractions and apply the Fourier inversion formula. We see then after some com- 

J CO 
e ilx j{x)dx, 
-OS 

where j(x) is given by 

(0, if z<0 






if x>0 



(2) 



3 National Bureau of Standards. 

2 University of Cincinnati. Deceased September 19, 1952. 

8 Figures in brackets indicate the literature references at the end cf this paper. 
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with 



C= 



aU(a 2 +bf)Ud 2 k 
j=i fc=i 



V U(a 2 + dl)Ubf 

k=l j=l 

771 / b 2 \ 



n). 



If we introduce the polynomial 



H(z)=Co(z-dl)(z-d 2 ) . . . (z-d£)=2>*z ro -*, 



where 



we can write 



<*>=(-!)■ 



m \ -1 

1=1 ) 



and c„ 



(2a) 

(3) 

(3a) 
(4) 



X 



\j=H(bf), j=l,2,...,n. 
The function /(«) as determined by (2) is real, after some computation it is seen that 

o 

f(x)dx—\. Iif(x) is nonnegative, it is the frequency function of an absolutely continuous 



probability distribution. We obtain therefore the following condition: 

Condition I. The rational function <j>(t) is a characteristic function if, and only if, 



9(e)- 



(1 — Xi cos bi6 



(1 — X 2 cos b 2 0), 



(I — \ n cos b n 6 



A2(n -1) 



7)2(7* -1) 



>0 



(5) 



for all 0. 

If we consider a rational function without zeros, that is, if we assume that O (O is the 
reciprocal of a polynomial whose roots are —ia, —iv j} —ivj (j=l,2, . . . , n), then we 
derive by a similar reasoning the following statement: 

Condition II. The function <£ (0 is a characteristic function if, and only if, 

1 — cos bid j 1 — cos b 2 6, . . . , 1 — cos b n 6 



goioy- 



«, 



bl 



7)2(71-1) 



7^2 (n -1) 



>0 



(6) 



for all 6. 

This condition is analogous to (5), in fact, (6) may be obtained from (5) by putting X ; =l 
for^l, 2, . . . , n. We agree, therefore, to put X^=l whenever m = 0. 

In this paper we consider the determinant (5) only for integer values of the bj and show 
that it represents for certain configurations of the bj and d k a nonnegative trigonometric poly- 
nomial. In this manner we can construct certain rational characteristic functions. 
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2. Some Lemmas 

For the discussion of the nonnegativity of g(ff) we need the following lemmas: 
Lemma 1. If b u b 2 , . . . , b n are integers, the determinant (5) can be factored so that 

0(0) = (1 — cos 6) n - m A(cos ff), (7) 

where A(x) is a polynomial in x of degree b n —n-{-m. 

Lemma 2. If b u b 2 , . . . , b n are odd integers, and if n^>m, we have 

g'(6)=sm 0(1 — cos 2 0) n " m - 1 5(cos 0), (8) 

where B(x) is a polynomial of degree b n —2n J \-2m J r \. 
Lemma 3. Let 

P s (x)= V {\-t 2 ) n ~ l t 2s dt fors=0, 1, (9) 

and 

Q(xS)=P 1 (x)-?Po(x). (10) 

The polynomial Q(x,£) is nonnegative in the interval — 1 <x<\ if, and only if, 

f<2n. (11) 

Here z n is the root of the equation 

R(z)^Q(-2,z)^j^(l-tY-\t 2 -z 2 )dt = 0, (12) 

which falls into the interval (0,1). 

Moreover, let w n == -J r 
Then lim w n =p exists and is the root of the equation 



n + 2 z ^ 



F(fi)=jT$-p^-2p*£{l-fy-**dt=0 



(13) 



which is located in the interval (0, Vi) • 

The proofs of the lemmas 1 and 2 are completely analogous to the proofs given for the 
corresponding lemmas in an earlier paper [3], in which nonnegative trigonometric polynomials 
of the form (6) were studied. Lemma 3 is identical with lemma 3 of [3]; we omit therefore 
the proofs of the three lemmas and refer instead to [3]. 

The function g(0) as defined by (5) and (2a) depends on the parameters 
bi, . . ., b n , di 9 . . ., d m . It is sometimes convenient to indicate this fact by writing 
g(6)=g(6\b 1 , . . ., b n ;di, . . ., d m ). Similarly, we shall write occasionally V =Vo(b h . . ., b n ) 
for the Vandermonde determinant formed from the numbers &?,...,&*. 

It is possible to obtain some of the results of our previous paper as particular cases of 
results obtained in the present paper. For this purpose we need the following lemma. 

Lemma 4. Let g{B\b x , . . ., b n )di, . . ., d m ) be the function defined by (5) and {2d). If 
we put d v =b p , we obtain 

g(0\bi, . . ., b n ) d\, . . ., dp-i, b p , d v +i, . . ., d m ) 
V (b u . . ., b n ) 

, Q\Q\b\, . . ., bp-\, 6p+i, . . ., b n ; d\, . . ., (z p -i, d v +i, . . ., d m ) , « * v 

Vo(6i, • • -, 6p-i, bp+i, • • •, K) 
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that is, the substitution d v = b p transforms the function 

g(6\b u . . ., b n ; d u . . ., d m ) 
V (6i, . . ., b n ) 

into a similar function that does not contain any more the parameters d v and b p . 

It is possible to give an elementary proof of (14) by a straightforward computation. 
The validity of (14) becomes apparent, however, if we remember that the rational function <j>(t) 
is the Fourier transform of Ce~ ax g(x\b h . . ., b n ; d u . . ., d m ). 



3. Nonnegative Trigonometric Polynomials 

In the following we assume that the b x , b 2 , . . ., b n are positive integers, while the 
di, . . ., d m may be any positive real numbers. The function g(6), denned by (5) is then 
a trigonometric polynomial. According to condition (I) the function <£(£) is a characteristic 
function if, and only if, g(6) is nonnegative. The purpose of this paper is to consider con- 
figurations that lead to nonnegative trigonometric polynomials. For this discussion it is 
convenient to introduce occasionally a new variable: 

x=cosd, (15) 

and to use the notation 

P(x)=g(&rc cos x). (16) 

Since we assume from now on that the 6/s are integers, the substitution (15) transforms 
cos bfi into T bj (x), where T b .(x) is the 6^th Tchebycheff polynomial of the first kind. There- 
fore, we see from (5) that under the conditions of lemma 1, 

In a similar manner, we have under the conditions of lemma 2, 

P\x)=ib (-l) j ^\ j T'(x) = -(l-x 2 y- m - 1 B(x). (18) 

pi J 

When we make this substitution we shall use the fact that the trigonometric polynomial g(S) 
is nonnegative for all if, and only if, the polynomial P(x) >0 for \x\ < 1. From (17) we obtain 
immediately 

Condition (III). The trigonometric polynomial g(0) is nonnegative if, and only if, the 
polynomial A(x), introduced by lemma 1, is nonnegative for \x\<l. 

We use condition (III) and the lemmas obtained in section 2 to derive conditions for certain 
trigonometric polynomials. The results are contained in the subsequent statements. 

Statement 1. If the b 3 are the first n consecutive integers, i. e., bj—j for j=l, 2, . . ., n r 
and if m=l, the trigonometric polynomial g{6) is nonnegative if, and only if, 2dl>n. 

In this case we have \ j = l—j 2 ld\, j=l, 2, . . ., n and obtain from (17) 

^(z) = g(- l)'- 1 ^ [l ~(l -^ 

From this equation it is possible to determine the coefficients a and 6. After some elementary- 
computations we obtain 
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Clearly, P(x)>0 for |z|<l if, and only if, \n 2 — d\\ <n(n— l)+d?. This inequality can be 
written as —n 2 +n—dt<n 2 —dl<n 2 —n J r d 2 . Since n is a positive integer, this condition reduces 
to 2d\>n. 

Corollary to statement 1. If the bj are the first n consecutive even integers, i. e., if 
bj=2j for j=l, 2, . . ., n, and if m=\, then the trigonometric polynomials g(ff) is nonnegative 
if, and only if, d\>_2n. The corollary follows almost immediately from statement 1. 

Statement 2. If the bj are the first n consecutive integers, i. e., if bj=j for j=l, 2, . . ., n, 
and if m = 2, the trigonometric polynomial g(6) is nonnegative, if, and only if, a certain polynomial 
Q(x) is nonnegative for |z|<l. Here 

Q(x)=L(n,d 1 ,d 2 )x 2 +M(n,d l ,d 2 )x+N(n,d 1 ,d 2 ), (19) 

where 

L(n,d u d 2 )=n 4 -n 2 (d 2 +dl)+d 2 d 2 2 , \ 

M(n,d 1 ,d 2 ) = 2n 4 -6n 3 +4n 2 -n+n(d 2 +d 2 )-2d 2 1 d 2 2 , (20) 

N( n ,d 1 ,d 2 )=n 4 -Qn' + 7n 2 -2n+(n 2 -n)(d 2 +d 2 )+d 2 d 2 2 . ) 

Proof. From (2a) we see that 

x ^( 1 -|)( 1 -|) 

for j=l, 2, . . . n, and obtain from (17) 

P (*)=g (- 1)^^ 1 1 -(l - jj) (l -|j) Tj(x)} = (l-xr-Xax*+bx+c). 

From this equation it is possible to determine the coefficients a, b, and c. After some elemen- 
tary computation, one obtains 

P(x)=~~^(l-xr-'Q(x). 

A particular case can be discussed easily. If rfi=(n+l) and d 2 =(n+2), we obtain from (20) 
i=4(n+l)(2n+l), M=-4:(2n+l)(2n 2 +n+2), N=4n 4 +4:n* + 19n 2 +5n+4. It is then easy 
to show that Q(x) is positive in — l<x<-\-l so that P(x) is nonnegative for the 
values kj=j(j=l, 2, . . . , n) and m=2, di=(n-\-\), d 2 =(n-\-2). This means in probabilistic 
language that the rational function with 2^+1 poles —i(a±ij)(j=l, . . ., n) and — ia and the 
four zeros — i(a±id{), —i(a±id 2 ) is a characteristic function. 

Statement 3. If the bj are the first n consecutive odd integers, i. e., if bj—2j—l 
for j=l, 2, . . . , n, and if m=l while ^>1 4 , then the trigonometric polynomial g(d) is non- 
negative if either one of the following two conditions is satisfied 



(i) d><y/2n—l, 



(ii) d<^2n-l, 
but 



v 



(2n-\)-d 2 
(2«-l) 2 -d 2 - "- 1 ' 



4 The case n = m=\ is covered by statement 1, and can therefore be excluded here. 
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where z n ~\ is the root of the equation 



1 (l-t 2 ) n - 2 (t 2 -z 2 )dt=0, 



which falls into the interval (0, 1). 



lj neither (i) nor (ii) is satisfied, then the junction g(Q) assumes also negative values. 
Coeollary to statement 3. If for large n and d<C2n—l 



4 



K) K2n-i)-^ p=0 _ 5939157 



(2n-l) 2 -d 2 
then g(6) is nonnegative. 

In formulating this statement, we wrote for d x simply d. From (2a) we see that 

X,= l-^^fori=l,2, 



n. 



According to lemma 2, as expressed by (18), we have 

P'(x)=^± (-i)^ 2 -(2j-i) 2 ]A,TVi(*)= -{i-x*y- 2 B{x), 

where B(x) is a polynomial of the second degree. From this equation it is possible to determine 
B(x), we obtain 



*-{'-¥H-^4 



(21) 



where 

T = 2 2w - 2 (2n-l)A rt >0. 
From (21) we see that 



B(-i)=B(i)- 2g "- ?,'"-'> ■,; 



am_[i-V] 



d? 

7>0 



(22) 



2^)=-2,[l-^]*; B>> {x) =-2y[l-^] 

We first assume d>2n—l, then B'(x) <0 for 0<x<l. Since 5(1) >0, we have i?(x)>0 in the 
interval (0,1), from B(— x)=B(x) it is seen that B(x)^>0 for |x|<l. Next we suppose 
(2n— \)^>d>-yJ2n— 1, then B(x) has a minimum at the point x=0, so that 5(x)>B(0)>0. 
Summarizing, we see that B(x)>0 for |#|<1, provided d>^2n—\. But this means that 
P'(x)<0 in \x\<l. From (17) it is seen that P(1) = if n>m, so that P(x)>0 in |x|<l, and 
therefore g(6) >0 for all 0, so that condition (i) is proved. 



We next consider the case d<CrJ2n—l. We write 



*K^-']4 



(2n-l)-c? 2 
~(2n-iy-d 2 
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If we denote by 

and by 

then we have 
Therefore, 



V-p-^-.>>o 



t 



_(2n-l)-d 2 



>0, 



~(2n-l) 2 -d 2 ' 
B{x)=y'{x 2 -?) and P'{x) = -y'{l-x 2 ) n ' 2 ix 2 -^ 2 ). 

P(x)=y' f\l-t 2 ) n - 2 (t 2 -y 2 )dt. 

Condition (ii) and the corollary to statement 3 follow then immediately from lemma 3 ; a brief 
table of the z n is contained [3]. 

We proceed to discuss a more complicated case by assuming that the rational function 
(/>(/,) has two zeros (i. e., m=2) and that it has the first n consecutive odd integers as poles, so 
that bj=2j— 1 tar j=l,2, . . . ,n. 

We obtain then from (18) 



where 



Wo write 5 



so that 



P'(x)= -7(l-x*)"-»[A+Ca , +Ite , ] > 
y=2 2 "- 6 (2n-l)A a 

C=4[(2n-9)\ n -(2n-l)\ H _ 1 ] 
A=(2n 2 -13n+22)\ n -(2n-l)(2n-7)\ n . l +(n-l)(2n-l)\„^J 



(23) 



(23a) 



G A 

4X re 



A 
"4X,' 



P'(x)=-4yK(l-x 2 ) n - 3 (4x i +ax 2 -b). 
For n^>m=2 we have P(1) = 0, and therefore 6 

P(x)=4y\ n r(l-t 2 )"-\4t*+at 2 -b)dt for m>3. 



(24) 
(25) 

(26) 



The nonnegativity of the polynomial P(x) depends on the nature of the roots of the polynomial 
B(x)=4:X' i -\-ax 2 — b. This introduces the need to distinguish a number of cases, for this purpose 
we write A(r)=4r 2 +ar— 6=4(r — ti)(t— t 2 ), where 



Tl =(-a-Va 2 +166)/8 
r 2 =(-a+^a 2 +Ub)/s 
B(x)=h(x 2 ) 

In the following we denote by 

R i (z)= P (\-t 2 ) n -H 2 «- l) (t 2 -z 2 )dt for {=1,2. 
We are now in a position to formulate statement 4. 



(27) 



5 We assume here Xn?*0. We need not consider the case X n =0, since this reduces hy lemma 4 to the case covered by statement 3. 
8 If n=*m, lemma 2 is not valid, and moreover, P(1)>0. The case n=m=2 would therefore require a separate discussion. 
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Statement 4. If the bj are the first n consecutive odd integers, that is, if bj=2j—l for 
j=l, 2, . . . , n, and if n^>m=2, then the trigonometric polynomial g(d) is nonnegative if, and 
only if, one of the following nine, mutually exclusive, conditions is satisfied 

(i) \ w >0, a>0, and b<0 

(ii) X w >0, a<0, and a 2 +166<0 

(iii) X„>0, a>0, 0<&<a+4, and R 2 (^t 2 ) — ti#i(Vt 2 ) >0 

(iv) X w >0, 0>a>— 4, 4+a>6>— a 2 /16, and B 2 ^t 2 ) — t 1 E 1 (^t 2 )>0) if 6<0, ^ Aaw *A« 

additional condition R 2 (— V r i)~~ t 2 Ri(— Vn) >0 

(v) X w >0>a+8>a+4>6>-a 2 /16 

(vi) X„>0, -4>a>-8, 0>a+4>&>-a 2 /16, and afeo R 2 (^t 2 )—t 1 R 1 (^t 2 )>0] if b^O, 

we have the additional condition R 2 {— ^t x ) — r 2 i?i(— V t i) >0 

(vii) X„<0, a>0, and 6>a+4 

(viii) X n <0, a<0, and b> max [0, a+4] 

(ix) X n <0, a+4<6<0, and B 2 (Vn) — T 2 Bi(VrO <0. 

The proof of this statement is rather tedious. It is, however, similar to the proof of state- 
ment 3 of the present paper and of statement (F) of [3] and is therefore omitted. 

We conclude this section by using lemma 4 to show that the earlier results may be obtained 
as particular cases of the statements of the present paper. We consider as an example 
statement 1 and write n=s J rl and put d x =k, where k is an integer. From lemma 4 we see 
that g{6\\, . . . ,s + l;k) and go(0\l, . . . ,k—l,k-\-l, . . . ,8+1) are simultaneously non- 
negative so that statement 1 reduces to statement C of [3]. 
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